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We have developed Bayesian formalism to describe the 
process of continuous measurement of entangled qubits. We 
start with the case of two qubits and then generalize it to an 
arbitrary number of qubits. 

The problem of quantum measurement of a qubit (two- 
level system) received renewed attention recently in re- 
lation to its importance for quantum computing. The 
case of sufficiently fast (instantaneous) measurement cau 
be readily described by "orthodox" collapse postulate,EJ 
and this is the case assumed at present by all quantum 
algorithms. However, in practice, especially for solid- 
state qubits, the act of measurement is not instantaneous. 
Because of typically low coupling between a solid-state 
qubit and a detector, it takes a considerable time be- 
fore the qubit state is completely destroyed by the act 
of measurement. Correspondingly, because of fundamen- 
tally unavoidable noise of the detector, the information 
about the state of measured qubit is available not imme- 
diately, but only after some time sufficient to get an ac- 
ceptably large signal-to-noise ratio. It is important that 
the timescale of measurement (and collapse) process may 
be comparable to the timescale of "free" qubit evolution 
(e.g., due to Rabi oscillations) or to the duration of the 
detector on-off operation sequence. (For example, if the 
detector is switched off when signal-to noise ratio is still 
on the order of unity, the measurement is only partially 
completed.) 

So, for practical needs we should be able to describe 
the measurement of a solid-state qubit as a continuous 
process. The formalism suitable to describe a continuous 
measurement of an ensemble of qubits has been devel- 
oped two decades agoB (for. its use in quantum comput- 
ing problems see, e.g., Ref.d). In contrast, the formalism 
describing the process of measuremejat of a single qubit 
have been presented only recentljuQ and is still in the 
stage of active development. (In fact, it can be consid- 
ered as a direct continuation of the well-developed field 
of selectivg ol conditional quantum rpjeasurements - see, 
e.g., RefsJjIlj and references in Ref.El). This formalism 
is called-JBayesian (because of essential role of the Bayes 
formulall3b and combines advantages of the "orthodox" 
approachEJ (ability to tgeat single quantum systems) and 
the Leggett's approachB (ability to treat continuous mea- 
surement). 

The Bayesian approach has been applied soJm'-qjiIjz. to 
the continuous measurement of a sing lequbit.HBlaliain 
this paper we apply it to derive the equations describing 
continuous measurement of entangled qubits. 
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FIG. 1. Schematic of two entangled qubits, one of which is 
continuously measured by a detector. The noisy detector out- 
put I it) is used to monitor the evolution of the double-qubit 
density matrix pit). 



Let us consider first the case of two entangled qubits, 
one of which is continuously measured by a detector (Fig. 
|l|) . As a main example, we consider qubits made of dou- 
ble quantum dots while detector is a quantum point con- 
tact (realizations based on single-electron transistors and 
SQUIDs are also possible - see Ref.El). Let us denote 4 
basis vectors characterizing the state of two qubits as 
|1) =|tT), |2) =|n), |3) =|iT), and |4) =\ii). (The 
basis for the first qubit is determined by its interac- 
tion with the detector, while for the second qubit the 
basis is arbitrary.) The qubits can interact with each 
other as well as be noninteracting (the entanglement 
can be a result of previous interaction). The free evo- 
lution of qubits is described by the Schrodinger equa- 
tion d'^/dt = {—i/fi)Hqb'^, where Hqi, is the Hamilto- 
nian of qubits only (not including interaction with the 
detector). Correspondingly, in the case without detector 
the density matrix p of double-qubit system evolves as 
dp/dt = {-i/D^lHqi,, p\. 

The detector output is characterized by two dc cur- 
rents, /| and 1 1, corresponding to two states of the first 
qubit, and the frequency- independent spectral density S 
of the detector noise. As usuaH we assume weakly re- 
sponding (linear) detector, |A/| <C lo, where A/ = 
and Jo = (J| -I- J|)/2, to neglect individual electrons pass- 
ing through the detector and consider the detector cur- 
rent I(t) as a continuous function of time. For the same 
purpose we assume that the timescale e//o (where e is the 
electron charge) is much shorter than other timescales in 
the problem (due to collapse, dephasing, and free evolu- 
tion of qubits). 

Let us start with the simplest case when qubits are 
"frozen", Hqi, = (so all the evolution is due to the 
measurement only), initial state of qubits is pure, and 
the detector is ideal (for example, quantum pioint con- 
tact at low temperature is an ideal detector,El as well 
e-electron transistor well inside the cotunneling 
). We can always represent initial pure state as 
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^ = a\^)(^ (at IT) + 6t li)) + /3 li) ® (a^ |T) + h ID): 
where the states of the second qubit are normaUzed, 
|a|p + = |a|p + — 1, and consequently 
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f3| = 1. Since the detector is not coupled to the 
second qubit, the evolution due to measurement affects 
only the factors a and wluch can be calculated using 
single-qubit Bayesian resulln (the overall wavefunction 
phase is of course not important): 
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where P| (r) and (r) characterize the conditional prob- 
abilities (for the first qubit in |t) and ||) states)-of getting 
a particular realization of the detector outputta I(t): 
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Pt(T) = {2TrD) -'-exp 
n(r) = (2vri?)-V2 
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In the language of double-qubit density matrix the evo- 
lution described by Eqs. (^ and (^) can be rewritten as 
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where p| = pii(O) -I- ^22(0) and pi = pssiO) + P44(0) 
correspond to initial probabilities to find the first qubit 
in It) and ||) states. 

If the initial state p(0) is not pure, its evolution can be 
calculated in the following way. Let us represent p(0) as 



p(0) = 5]p,(0)p,(0), 



(9) 



where ps is the classical probability of a pure state |s), 
Ps is its density matrix, and the sum is over a necessary 
number of pure states. (Of course, such representation is 
not unique in general.) To calculate p(r) we can apply 
"double Bayes" procedure: classical Bayes theorem to 
obtain probabilities Ps{t), 
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and the quantum Bayesian result [Eqs. (^-(||)] to cal- 
culate each Ps{t). It is easy to show that the resulting 



evolution of p(r) = ^gPs{T)ps{T) satisfies Eqs. (|1)-(||), 
which therefore are valid for arbitrary mixed states as 
well. Notice that Eq. (||) has an obvious interpretation 
as the conservation, of the "degree of purity" similar to 
the one-qubit casefl 

Besides the derivation of Eqs. (|1)~(||) using one-qubit 
Bayesian result (as above), they can also be obtained 
(in the case of pwcie states) directly using the "quantum 
Bayes tJaeorem" ,113 which says that the classical Bayes 
formulaic is applicable not only to the probabilities de- 
scribed by the diagonal matrix elements (that is obvious 
because of the correspondence principle), but also ap- 
plicable to the wavefunction amplitudes. Besides that, 
Eqs. (§)-(|l) can be easily derived "microscopically" in 
the case of a low-transparency quantum point contact at 
zero temperature. In this case, solving the Schrodinger 
equation for the|-qubitsj-jCoupled to the detector (for the 
model see Refs£2l andEl) one can obtain the following 
Bloch equationsEHI for the density matrix p^^ which con- 
tains index n corresponding to the number of electrons 
passed through the detector: 
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The equations for other components of are similar 
and can be obtained by the substitutions: p^^ P22, 
Ph P44, Ph P34> Pi3 P23' and ^5*3 p^^. Solving 
these equations and collapsing the number n at time r 
(measuring the charge passed through the detector and 
obtaining, for example, charge me), 



p7Ar + 0)=SnmP^J{T + 0), 
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one reproduces Eqs. (||)-(^). 

Now let us take into account finite detector j-jideal- 
ity (efficiency), rj < 1, where in one-qubit caseB rj = 
(A/)^4S'r is the ratio of the "information acquisition 
rate"E2l (A/)^/4S' and the ensepible dephasing rate F. 
Similar to the derivation of Ref.Q, let us consider first the 
case of a detector with neglected output (which is equiv- 
alent to "pure environment" ) . Then, averaging Eqs. 
(^ over the probabiHty distribution p|P|(t) -|- PxPx(t) 
of /(t) [see Eqs. (@)-(^)], we get the following: the right- 
hand side of Eqs. ^) and (^ becomes unity (which means 
that pii, P22, P12, P33, P44, and P34 do not change on av- 
erage), while the right-hand side of Eq. (^ is replaced 
by exp[— (A/)^/45'] (which means that pi3, pi4, p23, and 
P24 decay on average with the rate (A/)^/4S'). Similar 
to the one-qubit case, we can xegard a nonideal detector 
as two detectors "in parallel" ,E3 neglecting the output of 
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the second detector. In this way we obtain the following 
result for a nonideal detector: Eqs. (^) and (Q) remain 
valid, while Eq. (N) should be replaced by 



Piajr) ^ Pair) ^ P23(t) ^ P24(t) 

PuiO) PlM P23(0) P24(0) 



exp(-7r), 



(18) 



where 7 = (77-1 - 1)(A/)V45'. 

Notice that since the second qubit is not coupled to the 
detector, the state of the second qubit changes only due 
to its entanglement with the first qubit. In particular, 
in the case of no initial entanglement (when p(0) can 
be represented as a direct product), the state remains 
disentangled, the second qubit density matrix does not 
change, and Eqs. (^-(H) reduce to the Bayesian result 
for the first qubit. 

If the qubits are not frozen, Hgi, 7^ 0, the evolution 
due to Hqb should be added to the evolution due to 
measurement. Jn differential form [we use Stratonovich 
representation,^ so we take usual derivatives of Eqs. (|^), 
(^, and (p^)] we get the following Bayesian equations: 
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The equations for remaining components can be obtained 
from Eq. (|l|) by substitution {11} {22}, from Eq. 
(|^) by substitution {33} {44}, and from Eq. (H) by 
substitutions {13} {23} and {13} -> {24}. 

These equations allow us to monitor the evolution of 
the double-qubit density matrix if we know the initial 
state p(0) (for example, we have prepared qubits our- 
selves) and we know the detector output I(t) from an 
experiment. [To emphasize the noisy nature of I{t) we 
show this time dependence in Eqs. (p^"(|24|) explicitly, 
while the time dependence of the density matrix p is not 
shown explicitly.] To simulate the measurement process 
numerically, we need (similar to RefJj) to complement 
these equations by the formula 

Iit)-l0 = ^(pil+P22-p33-p44)+^(i), (25) 

where £^{t) is a zero-correlated ("white") random process 
with zero average and the same spectral density as the 




FIG. 2. Schematic of A'' entangled qubits, one of which is 
continuously measured by a detector. Bayesian formalism al- 
lows us to monitor A''-qubit density matrix p{t), using detector 
output I{t). 



detector noise, S*^ = S. [Eq. ( p5[ ) is derived from the 
probability distribution P'jP]'{T)+piPi(T) for the average 
current I(t) at sufficiently small r, so that evolution due 
to Hqb can be neglected.] 

The generalization of Eqs. (|l9|)-(p5|) to the case of ar- 
bitrary number N of entangled qubits, one of which is 
being continuously measured (Fig. ^), is pretty obvious. 
If both basis vectors i and j (from the set of 2^ basis vec- 
tors) correspond to the state |t) of the measured qubit, 
then the evolution of the matrix element pij is given by 
the equation 

Pt] = Y [-^9'" P^'^ + P-^ ~Y~ ^^^^^ ~ ' 

If both i and j correspond to the state | j) of the measured 
qubit, then 

Pi] = Y ^-^9'" P^'^ ~ - lo) ■ (27) 

Finally, if i corresponds to the state |t) while j corre- 
sponds to the state ||), then 

Pij = Y^Hqb,p\^j - (Pt - Pi)^{I{t)- Iq)- 7Py. 

(28) 

In these equations Hqb is again the Hamiltonian of qubits 
(without detector) while P|(t) and Pi{t) (now time- 
dependent) are the sums of the diagonal matrix elements 
of p{t), corresponding to the states |t) and ||) of the 
measured qubit. Eq. ( p5[ ) should be generalized as 

m - ^0 = ^ (/^T w - pm + m- (29) 

Now let us generalize the formalism to the case when 
the detector is coupled to all qubits (Fig. ||). Classi- 
cally, in this case there are up to 2^ different dc current 
levels li , corresponding to various combinations of qubit 
states. Some of these levels can coincide, for example, 
if the detector is not coupled to some qubits or if some 
qubits are coupled to the detector equally strong. Ap- 
plying the quantum Bayes theorem in the case of frozen 
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FIG. 3. N entangled qubits, continuously measured by a 
detector, coupled to all qubits. 

qubits, Hqb — 0, and taking into account finite ideality tj 
of tlie detector, we obtain the following equations: 



P,{t) = (27rD) -'"exp 



-1/2, 



Gxp(-7ijr), (30) 

-(l{T)^I,f /2d], (31) 
/4S', (32) 



where /(r) and D are defined by Eq. (||), and the sum 
in Eq. (BO) is over all 2^ basis vectors k (the basis is 
defined by the interaction between the detector and each 
qubit). Correspondingly, the probability distribution of 
/(r) is J2i PiiiO)Pi(l)- Notice that the exponent due to 
nonideality in Eq. (13^) disappears for diagonal matrix el- 
ements {i = j) and also if the classical currents li and 
Ij for two different configurations coincide. This is be- 
cause li = Ij means equal coupling of the detector to the 
states i and j, so the detector noise cannot destroy the 
coherence between these states. 

Let us briefly discuss what will happen to Eqs. (|30|)- 
( p^ ) if we relax the assumption of weak detector response, 

— Ij\ <^ {li + Ij)/2. As an example, let us again con- 
sider a low-transparency quantum point contact at zero 
temperature. In the case of moderate or strong response, 
each electron passed through the detector brings a sig- 
nificant information and, correspondingly, changes sig- 
nificantly the density matrix p of the qubits. Then the 
language of continuous detector current is not applica- 
ble anymore, and instead of considering average current 
/(r) we should count the number n of electrons passed 
through the detector during time r. Equation (3^) in 
this case does not change (except = since the de- 
tector is ideal), while the Gaussian distribution in Eq. 
( pl| ) should be replaced by the Poissonian distribution: 
PiXr) — (n!)^-'-(/ir/e)" exp(— /iT/e). It is not easy to 
introduce nonideality for a deteetor with finite response. 
If, however, we define 7y in a wayQ similar to optical quan- 
tum efficiency as a probability to observe an electron tun- 
neled through a detector (unfortunately, this definition is 
hardly justified in typical solid-state setups), then we can 
keep the exponential term in Eq. ( ^ ) and should replace 
Eq. (§) by = (^"^ - lKVT^ - /2e. 

Returning to the case of weak detector response and 
continuous current, differentiating Eq. ( |30| ) over time, 
and adding the free evolution due to Hq^, we finally ob- 
tain the following equation: 
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Equation (p5|) in this case is replaced by 
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Our final generalization is to the case of several detec- 
tors, coupled to N qubits. Each detector has its own set 
of up to 2^ classical current levels. It is important to 
notice that coupling of qubits to different detectors can 
define different sets of basis vectors. So, generalization of 
Eq. (^) requires to sum the terms due to measurement 
over all detectors, choosing particular basis for each de- 
tector. 

In conclusion, we have developed the Bayesian for- 
malism describing continuous measurement of entangled 
solid-state qubits. The case of two qubits, one of which is 
measured by a detector is considered in detail and then 
generalized to an arbitrary case. For nonideal detectors 
we have assumed the absence of correlation between out- 
put and backaction noises, so the formalism applicable to 
nonideal detectors with such correlationO still has to be 
developed. The results of this paper can be experimen- 
tally tested. However, such experiments seem to be still 
a little beyond the reach of the present-day solid-state 
technology. They could be attempted after proposed 
Bayesian experiments with a single solid-state qubit, in 
particular. Bell-type experiment.tj 
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